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Abstract 

We study Standard Model four-fermion production in e + e~ annihilation at LEP2 
energies and above using a semi-analytical approach. We derive the complete QED 
initial state corrections (ISR) to the reactions e + e~ — > (Z°Z°) — ► /1/1/2/2 and 

e + e- (W + W~) -► JififUi with h ¥= h and /< / e ± /z7 e ) . As compared to 
the well-known universal s-channel ISR, additional complexity arises due to non- 
universal, process-dependent ISR contributions from t- and u-channel fermion ex- 
changes. The full set of formulae needed to perform numerical calculations is given 
together with samples of numerical results. 



1 Introduction 



The LEP2 e + e _ accelerator will finally operate at energies between 176 and 
205 GeV [1] and thus pass the production thresholds for W ± and Z pairs. At 
LEP2, a typical process will be four-fermion production which is much more 
complex than fermion pair production as known from LEP1. This complexity 
shows up at tree level already, because of the many Feynman diagrams in- 
volved in four-fermion production. As tree-level amplitudes are not sufficient 
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to describe experimental data, radiative corrections to four-fermion produc- 
tion are needed. It would be desirable to derive complete 0(a) electroweak 
and QCD corrections to four-fermion production, but this has not yet been 
achieved, although many partial results were reported in [2-5] and in references 
quoted therein. In view of the anticipated experimental precision of LEP2 and 
future e + e~ colliders, it would be desirable to have theoretical predictions for 
cross-sections and distributions accurate at the level of half a percent. 

In recent years, three major approaches to four-fermion production in e + e~ an- 
nihilation have been developed, namely Monte Carlo approaches [3], the semi- 
analytical approach [6-13], and the "deterministic approach" [14]. Monte Carlo 
and deterministic techniques use numerical integration for all phase space vari- 
ables. Typically, the semi-analytical method performs analytical integrations 
over the five (seven, if ISR is included) angular degrees of freedom and uses 
high precision numerical integration for the remaining two (three, if ISR is 
included) squared invariant masses. It represents an approach to the high- 
dimensional, highly singular phase space integration problem inherent in this 
type of physical problem, which is elegant, fast, and numerically stable. Thus, 
it may serve as an ideal source of benchmarks for the other two approaches. 
From LEP1, we know that semi- analytical calculations are also relevant to 
experimentalists. 

In this article we will present semi-analytical results for the gauge boson pair 
production reactions 

CCS: e+e-^ (W + W~) ^f?f?ftf^), fi^M 
NC2: e+e"^ (Z°Z°) - fJ^Ul) , /i^/ 2 , fi^M 
NC8: eV-,(^° 1 A l7 7)-/i/i/2/ 2 (7) , /i^/a, frf^M (1.1) 

with complete initial state QED corrections. The three Feynman diagrams 
for the charged current CC3 process at tree level are given in figure 1. The 
two or eight diagrams for the neutral current NC2 and NC8 processes are de- 
picted in figure 20. Classifications of four-fermion processes may be found in 
references [8,15]. 

Initial state QED corrections (ISR) represent a dominant correction in e + e~ 
annihilation. Complete ISR to four-fermion production separates into a uni- 
versal, factorizing, process-independent contribution and a non- universal, non- 
factorizing, pro cess- dependent part. With the index J labelling the CC3, NC2, 
and NC8 processes, the ISR corrected cross-sections can be generically written 
as 



3 Strictly speaking, the NC2 process is well-defined (i.e. observable) only as on-shell 
reaction. 
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Fig. 1. The tree level Feynman diagrams for off-shell W pair production (the CC3 
process). Left: t- channel diagram. Right: s- channel diagram. The particle momenta 
are given by the hi and pj . 
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Fig. 2. XTie tree /ewe/ Feynman diagrams for off-shell neutral gauge boson pair pro- 
duction (NC8 process). Left: t-channel. Right: u-channel. The NC2 process is obtained 
by neglecting diagrams with exchange photons. 



da j. 



QEDl 



dsi ds 2 



J s 



(1.2) 



with invariant boson masses s\ and s 2 , reduced center of mass energy squared 
s', and tree level four-fermion production cross-section a Jt0 . The factor G(s'/s) 
contains all mass singularities Yn.{s/rri£) and incorporates the process-indepen- 
dent ISR radiators as known from s-channel e + e~ annihilation [16]. In addi- 
tion, there is a non-universal contribution o~j°Q^ mv . It appears together with 
t-channel and u-channel amplitudes and is mass singularity free. For pure s- 
channel contributions it is absent. 
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Complete initial state QED corrections were shortly communicated for the CC3 
process in [7] and for the NC2 and NC8 processes in [11]. While a definition of 
initial state radiation is straightforward in the neutral current process, there 
is an arbitrariness for W pair production in its definition. In [7], we restored 
the U(l)-invariance of the initial state photon emission by adding an auxiliary 
current. This arbitrariness is characteristic of charged current processes. Of 
course, it is the sum of the corrections what will be finally observable. 

In this paper, the complete analytical formulae for the non-universal correc- 
tions, supplemented by a study of their numerical importance, will be pre- 
sented for the first time. 

For on-shell production, 
e + e~ -> W + W-(-f), 

e + e" - Z°Z°( 7 ), (1.3) 
the generic cross-section may be obtained as follows: 



/ds' r 
— [G(s'/s)<tj, (s',M*,M*) 
s 



4M^ 



+ ^QErf V (s, s\My,My) . (1.4) 



Here, My represents the W or Z boson mass. The relation between aj : o and 
<7j ; o on one hand and and ^j°QED mv and (Tj°Q^ niv 011 the other hand will be 
discussed in section 2. 

For Z pair production, there are no additional QED corrections for the on-shell 
case, while for W pair production there are final state corrections and initial- 
final interferences. The non-universal initial state QED corrections were not 
known before as explicit analytical expressions. However, they have been deter- 
mined as a part of the complete electroweak corrections to the processes (1.3) 
with numerical integrations in [17]. 

The outline of this report is as follows. Section 2 presents general features of 
our approach to the complete initial state QED corrections. In section 3 we 
give a detailed presentation of the non-universal cross-section contributions. 
Section 4 contains numerical results and section 5 concluding remarks. In a 
series of appendices, we give technical details of the performed computations. 
Some notations are introduced in appendix A. Our phase space parametriza- 
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tion and all relevant relations between particle four-momenta and phase space 
variables are presented in appendix B. In appendix C, the tree level, the real 
ISR, and the virtual ISR matrix elements are given. These matrix elements 
represent the starting point for the calculations in this paper. The analytical 
integrals needed for the integration of the angular phase space variables and, 
in the case of virtual corrections, loop momenta are collected in appendix D. 



2 General Structure of Initial State QED Corrections 



To include initial state QED corrections (ISR) to (1.1), the five-particle phase 
space is required to take into account four final state fermions and a brems- 
strahlung photon with momentum p. We make use of the following parame- 
trization: 



1 V A ( S ' S '>°) V A ( s '' Sl ' S2 ) \/A(si,m?,mi) J X{s 2 , mj, ml) 

dl :, 



(2tt) 14 8s 8s' 8si 8s 2 

x ds' dsi ds 2 dcos^ dQ R dQ 1 dQ 2 . (2.1) 



In equation (2.1), the azimuth angle around the beam has already been inte- 
grated. We have adopted the usual definition of the A function, 



X(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2ac - 2bc, 

X = A(s, si, s 2 ). (2.2) 



We use k\ and k 2 as the initial electron and positron four-momenta, while pi, 
p 2 , ps, and P4 label the final state fermion momenta as indicated in figures 1 
and 2. The relevant invariant masses areQ 

s = -(k 1 + k 2 ) 2 = -(pi+p 2 +p 3 +p4 + p) 2 , 
s' = -(pi +p 2 +P3+ Pa) 2 , 
Si = -(pi +P2) 2 , 

s 2 = -(P3+P4) 2 . (2.3) 

The photon scattering angle 9 is defined as the angle between p and k 2 in the 

4 In our metric space-like four- vectors k have positive k 2 . Thus k 2 = —m 2 for 
on-shell particles of mass to. 
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Fig. 3. Generic two-boson production and decay Feynman diagram with real photon 
initial state radiation. 



center of mass system. The solid angle 

dQ R = dcos9 R d(f) R (2.4) 

represents the production solid angle of the boson three vector V1—P1+P2 in 
the four-fermion rest frame P1+P2+P3+PA — O. Oi [^2] is the solid angle of f>\ 
[p 3 ] in the two-fermion rest frame P1+P2 = [pz+Pa = 0]. In this frame, the 
three-vectors p\ and p 2 [P3 and p 4 ] are back to back. The polar and azimuthal 
decay angles in the above two-particle rest frames are defined via 

dtti = dcos#i d0i, % = 1,2. (2.5) 

Further details of the kinematics and the five-particle phase space may be 
found in appendix B. 

The processes (1.1) have the generic structure shown in figure 3. Their ISR- 
corrected cross-sections get two types of contributions, namely universal and 
non-universal ones. The first type is universal in the sense that it arises from 
photonic insertions to any of the basic diagrams and is independent of the de- 
tails of the subsequent interactions. These universal corrections are, of course, 
related to the collinear divergences of the (radiating) initial state electrons 
and positrons and appear to be known, including higher order terms, from, 
e.g., the study of ISR corrections to the (single) Z line shape. They are exactly 
those of the s-channel ISR contributions. For diagrams with t- and w-channel 
exchanges and interferences of these diagrams among themselves and with 
s-channel diagrams, there are additional corrections which are not logarithmi- 
cally enhanced and which depend on the details of the interfering amplitudes. 
These additional corrections are thus non- universal. 

ISR-corrected cross-sections for the processes (1.1) with soft photon exponen- 
tiation may be described by the ansatz 
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(2.6) 



with the threefold differential cross-section 



ds\ ds 2 ds' 



cf\ s ', Sl , S2 ) p^sy + ny 



Pe-l c( fe ) , <iy( fc ) 



(2.7) 



where 



2a 

7T 
f - 



In 



- f 



mi 



(2.8) 



In equation (2.7) we have used several additional notations which will now 
be explained. The subscript J labels different processes, J G {CC3, NC2, NC8} 
and the superscript index k stands for cross-section contributions which stem 
from squared amplitudes or interferences with distinct Feynman topologies 
and coupling structures. Using 



,(k), 



7TS Z 



(2.9) 



the soft+virtual contributions Sf^ and the hard contributions take the 
form 



Sj k \s, s'; si, s 2 ) = |l + 5(s)J af ,y)> (s'; si, s 2 ) + cr ( ?j(s'; si, s 2 ), 



H ( j\s, s'; si, s 2 ) = -H"(s, sQ cr { j k '°\s'; s u s 2 ) # + offi/s, s'; s u s 2 ) 



(2.10) 



Universal Part 



Non— universal Part 



with the well-known soft+virtual radiator S and the hard radiator H, 



fi( S ,. S ') = -i(l + fU + 0(a 2 ). (2.11) 



Higher order terms [16] may be implemented exactly as described in [2,3,13]. 
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Fig. 4. The amputated s-channel virtual initial state QED Feynman diagram. 





e e 
Fig. 5. The amputated s-channel initial state bremsstrahlung Feynman diagrams. 



We draw the attention to the definition of the (differential) effective tree level 
cross-section 



da Jfi (s'; si, s 2 ) y/X(s',si,s 2 ) 



dsi ds 2 



ITS 



12 



£cyv,*i,*2)0yv;*i,*2) (2.i2) 



which is inherent in the universal ISR correction part. Coupling constants 
and boson propagators are collected in cf \ while the Gp represent kine- 
matical functions obtained from the analytical angular integration. Explicit 
expressions for and Q^p may be found in appendix A. If the index k is 
associated with s-channel e + e~ annihilation only, non-universal ISR contribu- 
tions are absent. The s-channel ISR diagrams are generically shown in figures 4 
and 5. They contribute to W pair production. The non- universal cross-section 
contributions originate from the angular dependence of initial state t- and 
^-channel propagators and contribute both to W- and Z-pair production. 

To close the section, we comment on the on-shell cross-section (1.4). With the 
aid of equation (A. 14) it is obtained by the following replacements: 



a Jfi (s,My,My) = lim / dsi / ds 2 0j,o(s, si, s 2 ) 
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dsi / ds2 lim 

J T V ^Q 



S, Si, S 2 ) Gj \s, Si, S 2 



(2.13) 



and 



-non— umv / / 71 ^2 T\/r2\ 
^J.QED [S,S ,My,M v ) 

= dsi ds 2 V lim Cf\s', si, s 2 

J J , IV— +0 



non— umv / 



,(*) 



+ j(s,s';si,s 2 ) 



(2.14) 



3 Non-universal Initial State Corrections 



In this section, we present the final analytical results for the ISR corrected 
threefold differential cross-sectionsfor the processes (1.1). The starting point 
of our calculations are matrix elements of CC3 NC2 and NC8 processes, which 
are given in Appendix C. With the help of the computer algebra packages 
SCHOONSCHIP, FORM, and Mathematica [18], matrix elements were squared, 
spin summation was performed, the scalar products were expressed in terms 
of the phase space variables, and algebraic manipulations of the resulting 
expressions were carried out. All analytical integrations were obtained from 
hand-made tables of canonical integrals. The kinematical relations needed for 
the treatment of real bremsstrahlung may be found in appendix B. For the 
virtual corrections, tree level kinematics may be used as it was explained in 
appendix B of [13]. In the course of performing the various steps of analytical 
integrations we proceeded as follows. 

The virtual photonic corrections have been treated as a net sum of all con- 
tributing diagrams. The infrared singularity was isolated and, as a part of 
the universal corrections, subtracted from the net correction. Thus, the re- 
maining, non-universal virtual corrections are by construction free of infrared 
problems^. After tensor integration over the final state angular variables (see 
appendix D.l), the loop momentum integrations were performed with the aid 
of the integrals of appendix D.2 and the final integration over the vector boson 
production angle d in the center of mass system with the aid of appendix D.3. 

Also the real photonic corrections were first integrated over the final state 

5 The same was done with the real photonic corrections. Thus, non-universal virtual 
and real corrections may be treated completely separately. The interested reader 
may find details in [12]. 



9 




Fig. 6. The amputated t- and u-channel virtual initial state QED Feynman diagrams. 



angular variables (see again appendix D.l). Then we integrated over the pro- 
duction angles <pn and Or of the vector boson in the two-boson rest frame (for 
the canonical integrals see appendix D.4) and finally over the photon produc- 
tion angle 9 in the center of mass system (see the table of canonical integrals 
in appendix D.5). 

We have used the ultrarelativistic approximation for final state fermions and 
initial state electrons, i.e. the masses of these particles are neglected wherever 
possible. The various tables of canonical integrals were checked by Fortran 
programs with a precision of typically 10 -10 . At the end of our calculations, 
algebraic manipulations were carried out by hand to yield compactification of 
our final result S0 given in equations (3.3), (3.5), (3.8), (3.10), (3.14), (3.16). 

In the following two subsections, the non-universal corrections introduced in 
equation (2.10) will be explicitly given and commented. 

3.1 The Neutral-Current Case 

In the neutral current case, initial state QED corrections are represented by 
the Feynman diagrams of figures 6 and 7. 

The non-universal NC8 corrections are 

/ \ ^ ^1^2 non— univ ( \ 

(I , Oi S\S2 non— univ/ f \ /n i\ 

a Huc8( s ' s ^1^2) = 0fl, NC8 (s,s;sx,s 2 ), (3.1) 

7T 7TS 

6 The latter was cross-checked against FORM outputs with the aid of auxiliary FORM 
codes. 
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lip) 




Fig. 7. The amputated t- and u- channel initial state bremsstrahlung Feynman dia- 
grams. 



where the subindex V stands for virtual and the subindex R for real photonic 
corrections. 

The virtual contribution is a sum of contributions due to the different inter- 
ferences of loop diagrams with tree level graphs: 



non-univ 



°V,ncs {s;si,s 2 ) = <T vt 



non— umv 



, 11 non-univ . _non— umv 
+ 2a V,tu + a V,u 



The squared t- and w-channel contributions are equal: 



(3.2) 



non— umv 
T V,t 



(s; Sl ,s 2 ) = cr" iv (s; Sl ,s 2 ) 



3£ 2 
2 

+ A 



s L + 5 /_ — s (s — a) I 



12q 



, N . . 3 aS s — 4o" 

(s - 2a)(s - a) I 12q + -^j- Z_ — 1+ + s - a 

- ( A + a 2 + y J Il 2q - Y /_ - 4(7 /+ + 8(7 - 9s 



+ 2 [4 (2 - Z+) + (3s - (7) I 12q 



(3.3) 



Here, the following notations have been used: 



(7 = Si + S 2 , 
5 = Si - S 2 , 
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u = 



In^ + ln^ 



In— — In— = In— , 



In 



s s 

s — a + y/\ 
s — a — a/X ' 



Il9n — 



s (s — er) 
A 

g (s - o-) 3 
A 2 



Co - 



S2 



2V\ 



s — a 



2a/A 2 f V\ 



s — a 3 \ s — er 



4 (v + - i<_£. 



u " -A V 

S + 5 + a/A 



£12 = In 
£ 34 = In 



£>+ = Li 2 



s + <5- a/A ' 
s - 5 + \f\ 
s-5-y/X ' 

trnr 



Sl 



Li 5 



si 



+ Li 2 



S2 



Li- 



«2 



imm = - (s - er - a/a) , 

tmax = 7: (s ~ CT + a/A) • 



(3.4) 



The virtual corrections due to the tw-interferences are: 



_non— univ/ \ 
°V,tu (S, s l, s 2) 



3£ 2 



s Z+ + 5 1- — s(s — er) I 



12q 



+ C 



-a) ^l + 5sI 12q + ^/_^ - ^(sl + + 5lJ) 

7 -Cl-lQV d + AV d+ -A(l a -2l + + ^l a 



s — a 



+ 4Z_Z d _-4(l-^) - 3 (Z+ - 3) Z+ 



12 



+ 2 s 2 



(3 h 1+ + h 2+ ) (l a - -Z+) + - ( 3 In- + h 2 -) I. 



+ 2h 



1+ 



+ s 



9/ Rl 15 T 1 1 Q 

2 /+ - 6 l a - — s Iia, - 3- - 3- - 3 
2 di d 2 



+ a(l + -2l a ) + S -l\ 



+ 



8 s 2 



s — a 



+ s 2 (3 + h 2+ ) (V\I 12(l -2V a ) 

+ s y/X (2 I 12q + /n_Z_) + 2 (3s + a) V a 



(3.5) 



In equation (3.5), the following additional notations have been used: 



hi± 
h 2 ± 

d 2 
V d 

v d+ 



S — S\ s — s 2 



s 2 



s- s 2 

5 

Si 
S — Si 

J 

s 2 

s — a 



Sx 



(s - s 2 ) 



2 ' 



= In 



lndx — lnd 2 , 

tra.aa 



= Li 



Re 



s — a 



Li- 



t, 



S — (7 



LL 



Li 2 (di) + Li 2 (d 2 ) 



= L15 



s - s 2 

tmax 
S — S 2 



L19 



Li; 



train 



S - S 2 

train 
S - S 2 



+ Li 2 



Li; 



t, 



S — Sx 

traax 
S — S\ 



Li; 



+ Li 



train \ 
S — Sx/ ' 

train \ 

2\ , 

S — Sx/ 
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Li 3 
Li 3 



t r . 



Li, - 



s — a 

traax 
S — S 2 

trn.n 



— Lie 



LL 



tr. 



s — a 

train 
S - S 2 



+ Li 3 



S — Si 



Li, 



S — Si 



^1 tmir, 



- Li, - 



trnir 



+ Li 2 



d 2 i 



2 t-min 



- Li 3 



train 



dot 



2 L max 



(3.6) 



The non-universal real bremsstrahlung contribution to the NC8 process is given 
by: 



_non— univ / /. \ 
^RHCS l S ' S,Si,S2) 



_non— univ , _non— univ , _non— univ 



(3.7) 



Again, the squared t- and w-channel contributions are equal: 



non— univ 



'R,t 



1 a 



1 - 



,/2 



+ 



a (3s — s'_ — o^j 



A 



s a \ s si s 2 c 

H -1 + ^+12 



1 

+ 4 



A 



(7s' + S — (T 2s(7 2(7 ^ , 

Si (7 



A 2 



(L c i — L c2 ) 



ss' si 



:sa zct / 



+ 24 



A 2 



(L c i + L c2 ) 



7 + 



a 3asiS 2 



as's'_ 



As' 2 2A A 2 



(L c3 + L c4 ) + 2gg/2 £ c5 



4V 7 ! 



2(7 + s + si + j(^s [s'_ + a) - 4sis 2 ^ 

S Sl S 2 (7 (si + (7 



+ 12- 



A 2 



12 



14 



(a - s') s' ( , A 
+ J + l Dt i+ D l) + 



as', — s's'_ 



4s' 2 



D z \t2 + D z 2tl 



(3i 



In the above equation we have introduced the symbols: 



s' ± = s±s', X' = X(s',s 1 ,s 2 ), 



A = \(s'_, — si, — s 2 ) = si 2 + 2sV + 5 2 



(3.9) 



The additional notations are introduced in appendix D.5. For easy reference, 
the number of each integral corresponding to the enumeration in appendix D.5 
is listed in table 1. 

The most cumbersome contribution is due to the t-u-interference: 



non— umv^ I, 
T R,tu 



2 s' a 



s,s'; Sl ,s 2 ) = - {s -_ - } L c5 + -±(d\ + D\) 



s 2 + s' 2 



s 2 + s' 2 



~ 2sL(s'-a) ( D *™ + D *™) ~ (7 + 2sL(s'-a) J { D ™ +D '™) 
s 2 + s' 2 - Aa (s' - a) 



2 (s' - a) (4 - 2a) 



(L c6 + L c7 ) L c8 + D% + D 



s — 2s\ 



D 



'SSi 



t\u2 



s - 2s 2 

\/sS~2 



Dt2ul + 



s 2 + s' 2 



A sL (s' - a) 



+ 



+ 



s 2 + s' 2 - Aa{s' - a) s(s' -a)- Sl (s[ - 2a) 



2{s' + -2a)(s'-a) 
VX Is 2 + s' 2 - Aa{s' - a)] \s(s' -a)- s 2 « - 2a) 



2(4 -2a)(s'-a) 



AV- (3-10) 



Again, additional notations are introduced and the number of each integral 
corresponding to the enumeration in appendix D.5 is listed in table 2. 



3.2 The Charged- Current Case 



It is in order to mention that, for the CC3 process, the separation of initial and 
final state radiation is not unique. Since, in the t-channel contribution, there 
is electric charge flow from the initial state to the final state, electric current 
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Table 1 

Notations for o"^°°~ umv (s, s'; s±, S2). The numbers refer to appendix D.5 



Notation 


L c i 


L C 2 


La 


L C A 


L C 5 


D{ 


-^2 


D\ 2 


-Dzlt2 


D 2 2tl 


Integral No 


7) 


8) 


7) 


8) 


7) 


14) 


15) 


31) 


11) 


12) 



Table 2 

Notations for o-^~ umy (s, s'; s±, $2)- The numbers refer to appendix D.5 



Notation 


L C 6 


L C 7 


L C 8 


Dziti 


D z 2t2 




T) tu 

u a2 


Integral No 


18) 


19) 


18) 


10) 


13) 


18) 


19) 


Notation 


Alu2 


Dt2ul 




u t2u\ 


u t\u2 


u t2u\ 




Integral No 


36) 


37) 


38) 


39) 


40) 


41) 





conservation is violated and one faces the problem to find a gauge-invariant 
definition of initial state radiation. As a solution we proposed in [7] what we 
call the current splitting technique (CST). In brief, the CST splits the electrically 
neutral t-channel neutrino flow into two oppositely flowing electric charges +1 
and -1. Charge -1 is assigned to the initial state, charge +1 to the final state. 
This enables a gauge invariant definition of ISR with photon emission and 
absorption from the t-channel exchange particle. An auxiliary current is added 
to the naive charged current for this purpose in appendix C. Thus, the CC3 
t-channel receives the same ISR corrections as the neutral current t-channel. 
Of course, when performing a complete calculation, the final state corrections 
have to take into account opposite auxiliary terms so that the net auxiliary 
effect will vanish. 

The non-universal corrections are due to the interferences of the s- and t- 
channel diagrams of figures 4 and 6 with the corresponding Born diagrams 
and of the diagrams of figures 5 and 7 among themselves. These interferences 
show less symmetry than the interferences of the t- and w-channel diagrams 
of the neutral current case. However, as mentioned earlier, the pure s-channel 
corrections have no non-universal parts: 

4,cc3 = <4,cc3 = °- ( 3 - n ) 

Further, for the pure t-channel non-universal contributions to the CC3 process 
one finds, of course, expressions equal to those for the pure t-channel non- 
universal contributions to the NC8 process, namely 
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5,CC3 



\ ^ ^1^2 non— univ/ \ 

S',S 1 ,S 2 ) — — — o~ a V,t { s i s l, s 2), 

71 O 71 S 



"W*. ^ *o = l ^7 ^r univ ( s - ^ 



(3.12) 



with o-^ n " univ from equation (3.3) and <^ n ~ univ from equation (3.8). 

What remains to be calculated in addition are the virtual and real ^-interfe- 
rence contributions. 

The virtual correction is given by 



a S1S2 



71 8 7TS 2 



non— univ . 
T V,st 1 



s; si, S2) 



(3.13) 



with 



non— umv / 
T V,st \ 



s; si, s 2 ) = - V\ 7s + 3a + 5 I- + (Asa + 2A - 5 2 ) I 



12q 



+ A 



3a 2 75 2 



5sa + 3A - — — + \2sa + A - y I /+ + 5 + - ) L 



a 



(3.14) 



The real bremsstrahlung contribution reads 

a S1S2 



<7 S,cc3( s ' s '; Si ' S2 ) = - 



7T 8 7T 



„non— umv/_ _/. _ _ \ 



(3.15) 
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- 2 




(3.16) 



3.3 Discussion 

To close the section we will now discuss some formal features of our results. 
Non-universal cross-section contributions are analytically rather involved and 
contain many di- and trilogarithms. Although expected, it is noteworthy that 
interferences are more involved than matrix element squared contributions, 
which is especially true for the tw-interference. Looking at the integrals col- 
lected in Appendix D, one sees that some of them are very complicated (in 
appendix D.5 e.g. the integral 31) for the real t-channel correction and inte- 
grals 36), 38), 40) for the real tw-interference). However, the final answers for 
t-, u-, and st-contributions, both virtual and real, are remarkably compact. 
The virtual contributions contain only one true dilogarithm, which arises from 
the three-point scalar integral Ii2 q (see equation (3.4) and integral 15) in ap- 
pendix D.2). Both these virtual and real contributions do not contain true 
trilogarithms. In contrast, the ^-interferences are rather cumbersome. This 
may be traced back to the angular integrations involving products of different 
space-like fermionic propagators. Virtual t-u-mterferences contain several true 
trilogarithms and real bremsstrahlung tw-interference contributions exhibit 
complicated complex-valued dilogarithms. The latter, however, could indicate 
that some simplifications were overlooked. 

We mention that the four non-universal cross-section contributions to the CC3 
process became much more compact after inclusion of the auxiliary terms than 
without them. 

The attentive reader will have noticed that the virtual+soft non-universal con- 
tributions are evaluated at s' rather than at s. This is to avoid an unphysical, 
5-distribut ion-like concentration of non-universal virtual+soft ISR corrections 
at zero radiative energy loss. 

In equations (3.1), (3.12), (3.13), and (3.15), both reasonable threshold and 
high-energy properties may be observed. For \/~\ or v^A 7 approaching zero, 
non-universal corrections vanish, which may be verified by inspection of the 
explicit expressions. Thus also the ISR corrected cross-section vanishes for the 
kinematical zero of the tree level or the universally corrected cross-section. 
At high energies, the screening property may be verified to hold. By screening 
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we mean the assurance that cross-sections fall sufficiently fast with rising s. 
For the CC3 case, this is ensured in Born approximation due to an interplay 
between the s- and t-channel matrix elements and relies a certain relation be- 
tween the gauge couplings (see e.g. [19]). In the NC2 and NC8 cases, the same is 
achieved by an interplay between the t- and w-channel amplitudes leading to 
a cross-section proportional to the screening factor sis 2 /s 2 (the couplings of 
the two channels are equal) [8]. From the universal QED corrections, no addi- 
tional problems arise, because the interplay between the various contributing 
tree level matrix elements is not disturbed. This is completely different for the 
non-universal corrections. Here, the various interferences get different correc- 
tions and the interplay between them is destroyed. Fortunately, one may see 
from the explicit expressions that the screening factor observed in the net NC2 
and NC8 tree level cross-sections is rediscovered not only for full cross-sections, 
but even for individual non-universal contributions. So, unitarity is not vio- 
lated. We were not able to find the screening property for non-universal CC3 
corrections without including the auxiliary corrections. Also in this respect, 
auxiliary corrections seem to be a quite natural, if not necessary ingredient of 
the calculation. 



4 Numerical results 



To illustrate the analytical results of sections 2 and 3, we present numerical 
results [20]. The effect of ISR on the boson pair production processes NC2 
and CC3 is seen from figures 8 and 9 where inclusive cross-sections are pre- 
sented. From both figures one realizes that, as discussed in section 2, the 
dominant part of the initial state QED corrections originates from universal 
ISR. Universal ISR is typically of 0(10%). Non-universal ISR on the contrary 
is suppressed and rises from a few parts per thousand in the LEP2 energy 
region to 0(1%) at 1 TeV. For the NC2 process at 1 TeV the non-universal 
contribution to the cross-section is 2.5%, whereas for CC3 at 1 TeV it is 1.4%. 
For reference we have given cross-section values for both Z and W pair pro- 
duction in table 3. The numerical precision is better than 10~ 4 . As numerical 
input for table 3 we have used the standard LEP2 input (table 5 in [3]) which 
is reproduced in table 4. The weak mixing angle is determined by 





(4.1) 



y/2G F M% 



and we use the relation GjV2 = g 2 /(8M^). 
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200 300 400 500 600 700 800 900 



Vs (GeV) 

Fig. 8. The inclusive total off-shell Z pair production cross-section af a f(s). In the in- 
set, the relative deviations of the universally and completely ISR corrected cross-sec- 
tions from the tree level cross-section are given. 

In figure 10, we present, as an example for the NC8 process, the cross-section 
for 

e+e - _» {Z °Z°, z o 7) 77) ^ ^- hl (7) _ (42) 

Due to the negative slope of the cross-section curve ISR corrections are al- 
ways positive. Universal ISR corrections vary between approximately 12% at 
^=130 GeV and 21% at 600 GeV. Non-universal corrections rise from 0.9% 
at 130 GeV to 4.2% at 600 GeV. Thus, compared to the NC2 case, the NC8 
non-universal corrections are considerably enhanced. 

Finally it is worth mentioning that non-universal ISR tends to be harder than 
universal ISR, a fact which emerges from a numerical analysis of cross-sec- 
tions with lower cuts on s' [12]. This may be understood as follows: 0(a) non- 
universal corrections are infrared finite, whereas the infrared divergent uni- 
versal corrections are resummed by means of the soft photon exponentiation. 
Thus the universal corrections contain important soft resummed parts, and 
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Fig. 9. The inclusive total off-shell W pair production cross-section cr^,^{s). In 
the inset, the relative deviations of the universally and completely ISR corrected 
cross-sections from the tree level cross-section are given. 



the non-universal corrections do not. 



5 Concluding Remarks 



In this paper we have presented the analytical results of the first complete 
gauge-invariant calculations of initial state QED corrections to off-shell vector 
boson pair production in e + e~ annihilation. We have decomposed the correc- 
tions into a dominant universal contribution and a numerically smaller non- 
universal contribution. Non-universal corrections were found to be smoothly 
behaving around the two boson production threshold and to be explicitly 
screened from unitarity violations at high energies. As a by-product, the an- 
gular distribution of radiated photons is available from intermediate steps of 
the calculation but has not been presented here. 

By construction, we have not taken into account genuine weak corrections [22] 
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Table 3 

Inclusive total off-shell cross-section values for the processes NC2 and CC3 with and 
without initial state QED corrections 





_tree 
W NC2 


ONC2 [pb] 

_umv 
°NC2 


°NC2 


-.tree 
W CC3 


OCC3 [pb] 

°CC3 


W CC3 


161 


0.0154 


0.0127 


0.0127 


4.8087 


3.4521 


3.4648 


175 


0.0648 


0.0506 


0.0508 


15.9168 


13.2903 


13.3408 


183 


0.3811 


0.2733 


0.2743 


17.6819 


15.4036 


15.4637 


192 


0.9783 


0.7788 


0.7817 


18.4479 


16.5958 


16.6625 


205 


1.1833 


1.0274 


1.0313 


18.5078 


17.1879 


17.2600 


500 


0.4097 


0.4424 


0.4479 


7.3731 


7.9168 


7.9810 


800 


0.2124 


0.2388 


0.2435 


3.9971 


4.4501 


4.5026 



Table 4 

Input parameters for table 3 



Quantity 


Value 


Quantity 


Value 


M z 


91.1888 GeV 


a,(0) 


1/137.0359895 


Tz 


2.4974 GeV 


a{2M w ) 


1/128.07 


M w 


80.23 GeV 




1.16639 • 10" 5 GeV" 2 


T w 




VcKM 


1 



or QED corrections related to the final state [17]. Compared to ISR, both are 
known to be smaller but may reach several percent nevertheless. Therefore 
they are comparable in size to the non-universal corrections calculated here. 
See also references [2,23]. Special suppressions of the order (9(1~V/My) are 
estimated for interference effects between initial and final state radiation in 
resonant boson pair production and inter-bosonic final state interferences [24]. 
The complete photonic corrections to four-fermion production are not available 
although important steps towards their calculation have been achieved [25]. 

The presented calculation should be considered in connection with the pro- 
gram of semi-analytical treatments of complete processes at tree level [8]. It 
would be quite interesting to compute complete ISR not only for the off-shell 
production of vector boson pairs, but also for true four-fermion processes. Sim- 
plest are the so-called NC24 [9] and CC11 [13] processes, but also the combi- 
nation of NC24 with Higgs production [10] should be straightforward. Because 
in these reactions all background Feynman diagrams, i.e. those diagrams that 
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Fig. 10. The total off-shell cross-section for process Jf..2. The inset shows how Zand 
photon pair production contribute to the NC8 cross-section. For the particle masses 
we used numerical values from [21]. 



must be added to CC3 [NC8] to get the sets for CC11 [NC24], are s-channel di- 
agrams, there is nothing but universal ISR for the pure background. So the 
only missing pieces are the non-universal ISR contributions to the interfer- 
ences between the boson pair signal and the background. The results may 
be of numerical relevance in the NC24 case below the ZZ production thresh- 
old. For CC11 and for the inclusion of Higgs production they will be of mere 
theoretical interest. 
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A Couplings and Factorizing Kinematical Functions 



A.l Couplings 



The electroweak couplings and neutral boson propagators enter the Born and 
QED corrected cross-sections in certain combinations. 

The NC8 process is described by one coupling function [9]: 



-NC8 



(6tt 2 ) 2 



V Re - 

v^tv^z D Vi {s,)D Vj (s 2 )D* Vk ( Sl )D* Vi (s 2 ) 

x [L(e,Vi) L(e,V k ) LicVj) LicV) 

+ R(e, Vi) R(e, V k ) R(e, V 3 ) R(e, V t ) 
x [L(/i,K> Vfc) 

+ R(fi,Vi) R(f u V k )] N c (f 1 )N p (V i ,V k ,m 1 ,s 1 ) 

x [LfaVj) L(f 2 ,V t ) 

+ RifaV^RihM)] N c {h)N p {Vj,V h Tn2,S2). 

{A.l) 



We have used the left- and right-handed fermion-boson couplings 



Hf, 7) = L(f, Z) = -J- (it - s* w Q f ) , 

R(f,l) = R(f,Z) = ~^-s 2 w Q f (A.2) 

z z Cw 

where Qf is the fermion charge in units of the positron charge e, // is the 
fermion's weak isospin third component, and sw denotes the sine of the weak 
mixing angle. Q e = — 1 and the fine structure constant is given by e = gsw = 
^/Ana . In addition, equation (A.l) contains boson propagators 

D v (s) — s — My + isT v /M v . (A.3) 

The photon has zero mass and width, M 7 = T 7 = 0. The color factor N c (f) is 
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unity for leptons and three for quarks. The phase space factors 

!1 for V m ^ 7 or 14^7 

(A.4) 
for v m = v n = 1 

in equation (A.l) correctly take into account that fermion masses may be 
neglected if the corresponding fermion pair couples to a Z resonance but may 
not if it only couples to a photon. 

The cross-section for the NC2 process is obtained by omitting exchange photons 
in equation (A.l). 

The CC3 process is described by three coupling functions 0$%, one for the 
i-channel, one for the s-channel, and one for the st-interference contribu- 
tion [6,13]: 



CC3 



(6tt2)2 \D w ( Sl )\ 2 \D w (s 2 )\ 2 

xL 4 (e, W) L\f u W) N c {h) L 2 (f 2 , W) N c (f 2 ), 



(A.5) 



C s 



CC3 



Cst 
i 



CC3 



(6vr 2 ) 2 Vi y j=l! z 



D w (s 1 )\ 2 \D w (s 2 )\ 2 D Vi (s)D* v .(s) 



X 



L(e, Vi)L(e, Vj) + R(e, VjR(e, Vj) 



x gsiVMVj) L 2 {h,W)N c {h) L 2 (f 2 ,W)N c (f 2 ), 

—^i Y Re 2 2 L 2 (e,W) 

(6vr 2 ) 2 V X Z \D w ( Sl )f \D w (s 2 )\ 2 D Vi (s) 

x L(e, Vi) gsiVi) L 2 (f u W) N c (f\) L 2 (f 2 , W) N c (f 2 ). 



(A.6) 



(A.7) 



The left-handed couplings L(f, W) of a weak isodoublet with the fermion / 
to a W boson are given by 

L(f.W) = ^ (A.8) 

We neglect Cabibbo-Kobayashi-Maskawa mixing effects. The right-handed 
couplings R(f, W) vanish. The couplings ^(V^) originate from the three-boson 
vertices and are 



5-3(7) = QSw, 
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g 3 (Z) = gc w . (A.9) 

Using G^/V2 = g 2 /(8M^) in the C functions, one easily rewrites the cross- 
sections for heavy boson pair production, namely the processes CC3 and NC2, 
in terms of Breit-Wigner densities 

»M-h«-wAw* Bm ' v=z ' w± (A ' 10) 

with branching fractions BR(i) into the appropriate fermion pairs. 
For the process NC2 one gets 

Pz(s 1 ) P z(s2) (G,M 2 ) 2 



CmC2 — 



(l-2s 2 w y+(2s 2 w " 4 



(A.ll) 



SiS 2 4 

For W ± pair production, the situation is slightly more involved. Defining 

(C M 2 \ 2 

C CC3 = W) Pw(si) P w(s 2 ), (A.12) 

sis 2 

one obtains by inspection of equations (A.5)-(A.7) 



^CC3 ~~ ^CC3 C uu , 

4 

Qc3 = ~CcC3 (c 77 + C^z + Czz) j 

s 

Qc3 = ~CcC3 (<Vy + Cvz) (A. 13) 

s 



with factors c ab as defined in [8] (see also [6]). 

For on-shell heavy boson pair production, T v —>0, Breit-Wigner densities are 
replaced by ^-distributions, 

Pv(s) ^ 5(s - My) x BR(i). (A. 14) 
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A. 2 Kinematical Functions for the Factorizing QED Corrections 



The kinematical function for the NC case is 



Qm2{ s 'i Sl, S 2 ) — ^NC8(s; Si, S 2 ) — S X S 2 



For the CC case, one finds [6]: 



s 2 + (si + s 2 f 
s — Si — s 2 



C-2 



(A.15) 



^CC3( S ' S l> S 2) 
^CC3( S ! S l? S 2) 



^cc3( s ; s i> s 2) - 



A 

192 



A + 12 s (si + s 2 ) -48sis 2 

+ 24 (s - si - s 2 ) S!S 2 C 

A + 12 (ssi + ss 2 + s x s 2 ) 

A + 12(ssi + ss 2 + sis 2 ) 
- 24 (ssi + ss 2 + sis 2 ) sis 2 C 



48 H* " *i " *2 



(A.16) 
(A.17) 

(A.18) 



for the kinematical functions. The logarithm C contained in integrated t- or 
M-channel contributions is defined by 



r rl v 1 s-s 1 -s 2 + \^X 

C = £{s;s 1 ,s 2 ) = -= In = 

V A s — Si — s 2 — V A 



(A.19) 



B The 2^5 Particle Phase Space 



When considering initial state photon bremsstrahlung to a 2^4 process, a 
photon with momentum p appears as fifth particle in the final state. The in- 
termediate vector bosons (two-fermion systems) have momenta V\ = pi + p 2 , 
V2 = P3 + P4- The five-particle phase space has eleven kinematical degrees 
of freedom. It is convenient to parametrize the five-particle Lorentz invariant 
phase space by a decomposition into subsequently decaying particles: 



1 \A(s,s',0) x/\{s',s 1 ,s 2 ) yJ\( Sl ,ml,ml) y/\(s 2 ,ml,ml) 
5 ~ (2tt) 14 8s 8^ 8s^ 8^ 

x els' dsi ds 2 dcos# d(j)R dcos#R dfa dcos^i d0 2 dcos# 2 , (B.l) 
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where nii is the mass of the final state fermion with momentum pi and 
: Azimuth angle around the beam direction. 

9 : Polar (scattering) angle of the photon with respect to the e + 

— * 

direction k 2 in the center of mass frame. 
(p R : Azimuth angle around the photon direction p. 
9r : Polar angle of v\ in the iyi+v 2 ) rest frame with z axis along p. 
<pi : Azimuth angle around the V\ direction. 

9\ : Decay polar angle of p\ in the V\ rest frame with axis along v 1 . 
02 : Azimuth angle around the v 2 direction. 

9 2 : Decay polar angle of ps in the v 2 rest frame with axis along v 2 . 

Si : Invariant mass of the final state fermion pair f±f 2 : s± = —vf. 

s 2 : Invariant mass of the final state fermion pair fsf4 ■ s 2 — — v^- 

s' : Reduced center of mass energy squared. This is the invariant 
mass of the final state four-fermion system: 
s' = -(vi + v 2 ) 2 = -(pi + p 2 + p 3 + p A ) 2 . 



The ranges of the kinematical variables are 




(mi + m 2 + m 3 + m 4 ) 2 


< 


s' 


< 


S 


(m 1 + m 2 f 


< 


Si 




V^s 7 — rn 3 — m 4 


(m 3 + m 4 ) 2 


< 


S2 


< 


(Vs 1 - v^r) 2 


-1 


< 


cos9 


< 


+1 





< 


0{fl,l,2} 


< 


2tt 


-1 


< cos6» {Rili2} 


< 


+1 



or, alternatively, 
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fl(pl) f,(pj) 

Fig. B.l. Graphical representation of a 2^5 particle reaction. In general, the vectors 
vi,V2,Pi,P2,P3, and pi do not lie in the plane of the picture. In the figure, the angle 
9 is called # 7 to emphasize that it belongs to the direction of the radiated photon. 

(m 1 +m 2 ) 2 < 
(m 3 + m 4 ) 2 < 

-1 < 

< <f>{R,l,2} < 2tt 

-1 < cos0 {jR;1 , 2} < +1 . (B.3) 

An illustration of the 2 — > 5 particle phase space is given in figure B.l. 

For the processes under consideration, the integrations over the final state 
fermion decay angles <f>i, 9 1 and <f> 2 , 6 2 are separated from the other inte- 
grations and are carried out analytically with the help of invariant tensor 
integration (see appendix D.l). The remaining integrations are nontrivial. We 
will now explain, how the scalar products appearing in the squared matrix ele- 
ments can be expressed in the center of mass frame with a Cartesian coordinate 
system as drawn in figure B.l using the phase space variables introduced in 
equation (B.l). 





< 1 


[y/s — m 3 — m A ) 


S2 


< 


(y/s- v^l) 2 


s' 


< 


s 


COS0 


< 


+ 1 
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The initial state vectors k\ and k 2 are given by 



ki = (k°, k sin#, 0, — k cos#) , 
k 2 = (k°, — k sin#, 0, k cos#) 



(B.4) 



with 



k° 



2 ' 



(B.5) 



where m e is the electron mass. 

For the momentum of the photon radiated along the z axis one finds in the 
center of mass system: 



p = (p°, 0, 0, p°) with 



P° = 



s — s 



(B.6) 



2v^ ' 

Thus, the invariant products k±k 2 and k(p may be expressed easily. 
The virtual boson pair recoils from the photon in the center of mass system: 
fs + s' 



(Vl + v 2 ) 



V 2 v^ 



, 0, 0, - 



s — s 



B.7) 



The vectors v\ and v 2 in the two-boson rest frame (the i?-system with v^r+ 
v 2 ,r = 0) are easily derived: 



vi,r= (pff,0AV R ) 



^S' + Sl -S 2 y/X(s',S 1 ,S 2 ) S ' 

, u, u, 



V 



2y/? 



2y/? 



v 2 



R = (p§ o ,0,0,-P fl ) 



/ 



S' - S 1 + S 2 \jK S 'i S l> S 2) 

, u, u, 



V 



2V? 



2V? 



. (B.8) 



Both vi t R and v 2 ^ R do not depend on angular variables. 



We need Vi in the center of mass system. The boost between the center of mass 
frame and the two-boson rest frame R with zero component 7^ and modulus 
of the spatial components 7 is 



30 



/rx\ S ~\~ S S S 



Taking into account that the production angles of the virtual vector bosons 
are not fixed, equation (B.8) yields the following center of mass system four- 
vectors: 



vi = (l {0) Pi2° ~ lV R cos6 r , V r sin6 R cos0 R , - V R sin6 R sin(f) R , 

^V R co S 9 R - 1P ?f), 

v 2 = ( 7 (°) p**f + 1 <p R COS R , - V R sin9 R cos0 R , V R sm9 R sin(/) R , 

-l i0) V R cos6 R - lp Kf). 

(B.10) 



The sum (vi jR + V2,r) in the i?-system is simply (a/s', 0, 0, 0). The boost into 
the center of mass system yields (v\ + V2) = (7 < - - ) \/s', 0, 0, — 7\/s'), which is 
exactly (B.7). 

The above formulae allow to express the products k(Vj and pvj. The product 
P1P2 [and analogously P3P4] is most easily calculated in the center of mass 
system i?x [R 2 ] defined by P1+P2 = [P3+P4 = 0]. This is explained in equa- 
tions (B.5)-(B.ll) of [13] for the 2—^4 process. The discussion of reference [13] 
remains valid also with an additional photon in the initial state. For the re- 
maining scalar products ppi, kjPi,VjPi, and piPj we used explicit expressions 
for the vectors pi in the center of mass system. The vectors pi are obtained 
from the vectors pi tR in the i?-system by rotating and applying the Lorentz 
boost (B.9): 



Pi 



7 



pIr( s ') cos °r - pIr( s ' 



sin9 



'R 



cos(p R + pl R (s') sin(p R 



i {0) pU s 

pIrW) wd9 r + pIr{s') cos0h 
PlA 3 ') sm9 R + pl R (s') cos9 r ] smfa + pj fl (s') C os(p R 
pIrW) cosd R - pIrW) sin ^] - IPiA 3 ') J 



V 7 



,(0) 



fB.ll^ 



The four-momenta PiA s ') are exactly the final state momenta pi shown in 
equation (B.21) of [13], when these are evaluated at s' instead of s. The frame 
R there has to be understood as frame R\ or R 2 introduced here. 



Now, the derivation of all momenta in the center of mass system in terms of the 
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integration variables is completed for the 2— »5 process and all scalar products 
appearing in the initial state bremsstrahlung calculation can be expressed in 
terms of phase space variables. 

For the calculation of virtual and soft photon corrections, the 2— >4 phase space 
is needed together with the corresponding momentum representations. It may 
be taken over completely from appendix B of [13]. 



C Matrix Elements 



In this appendix, we will present the matrix elements for the CC3 and NC2 
processes. The matrix elements for the NC8 process are easily obtained from the 
NC2 matrix elements by adding amplitudes where Z couplings (propagators) 
have been replaced by photon couplings (propagators). 



C. 1 Tree Level Amplitudes 



For the calculation of soft and virtual corrections, the tree level matrix ele- 
ments for the CC3 and NC2 processes are given in terms of s-channel, t-channel, 
and li-channel amplitudes corresponding to the Feynman diagrams of figures 1 
and 2: 



MiL = M* Z + M?^ + M B 



l CC3 
l NC2 



AC, = Mf z + M 



8,1 

B A aB 

u,Z 



t,Wi 



(C.l) 
(C.2) 



with 



M 

Mf w 



9n' 



9(3(3' 



9o 



Dz(s) D w (si) D w (s 2 



'-T Kz M?2, w M^ w g 3 (Z) (C.3) 



#77' 



9m 1 



9c 



DJs) D w (si) D w (s 2 



~T B]' M?' w M* gM (C.4) 



9(313' 



9o 



1 



Dw(si) 


D w (s 2 ) 


9/3(3' 


9aa' 


D z ( Sl ) 


D z (s 2 ) 


9aa' 


9(3(3' 



a(3 A ,/3' 



CM: 



M a 



D z ( Sl ) D z (s 2 ) ql 



't,Z 1V1 12,Z lvl 34,Zi 



'u,Z 1V1 12,Z 1V1 34,Z- 



(C.5) 
(C.6) 
(C.7) 
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We have used the vector boson propagator denominators Dy defined in equa- 
tion (A. 3), the trilinear coupling g^iV ) given in equation (A. 9), and 



(C.8) 



In the massless approximation, the annihilation matrix elements B] y0 , B" v , 



B^y, and the decay matrix elements M^ v are given by: 



B]' y0 = u(-k 2 )Y [L(e, V°) (1 + 75) + R{e, V°) (1 - 7 s)] u(h), (C.9) 

B# = 2 ui-h)^ [L 2 (e, V) (1 + 75) + R\e, V) (1 - 75)] it ^u(h), 

(CIO) 



B^ z = 2u(-k 2 ) 1 a L 2 (e,Z)(l + l5 ) + R 2 (e,Z)(l- l5 ) 4 u ^u(h), 



M*y = u(p t )r [L(fi, V) (1 + 75) + R(f» V) (1 - 75) 



(C.11) 

u(- Pj ) (C.12) 



with left- and right-handed couplings as defined in equations (A. 2) and (A. 8). 
We have further used 



q t = h-vi = -k 2 + v 2 = - (h - k 2 - V! + v 2 ) , 



q u = h-v 2 = -k 2 + vt = -(h - k 2 + v ± - v 2 ) 



(C.13) 
(C.14) 



C.2 Virtual Initial State Corrections 



The amplitudes for virtual initial state corrections are easily obtained from 
the matrix elements in equations (C.3) to (C.7) by substituting the initial 
state currents. 

The t-channel NC2 virtual QED matrix element Aif z is represented by the 
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Feynman diagrams of figure 6 and given by 
M 



V _ 9/3/3' 9aa' T/a/ 3 , W3' W 

l ' Z ~ D z ( Sl ) D z (s 2 ) V t,z M i2,zM u , z 



(C.15) 



with 



| 16tt 2 



- i // 4 -«) 



/■ d n p / 
7 (2tt)» V 



v t,vertl + 17 * 



, , ya/3 \ 

t,vevt2 ^ v t,setf ^ v t,box J 



(C.16) 



fj,=m e 



where dimensional regularization is applied with p being the loop photon mo- 
mentum, and where 



ya(3 



t,vertl 



v a/3 
v t,vert2 



u(-k 2 ) 



-2^-7^)7° (it-f)lJtlP 



p 2 — ie 




(P - <?t) 2 - i£ 




p 2 + 2/c 2 p — is 




x2 


\L 2 (e,Z) (l + l5 ) + R 2 (e,Z) (f 


-7s) 



u(-k 2 ) 



Tin 11 {it-i>)i p {2k^-i> llx ) 



p 2 — ie 




(P - (It) 2 - ie 




p 2 — 2k\p — ie 


ql 




x2 


\L 2 (e,Z) (l + l5 ) + R 2 (e,Z) (f 


-75) 


u(fci) 



t,self 



u(-k 2 ) 



Titi 11 (4t- 



p 2 — ie 

x2 



(P-Qt) 2 - ie] [q 2 
L 2 (e,Z) (l + l5 ) + R 2 ( e ,Z) (1 - 75 )1 u^) , 



K t,box 



u(-k 2 ) 



x 



:-2A#- 7 "j07 a {it~i)l P (2^- 



p 2 — ie p 2 — 1k\p — ie p 2 + 2/c 2 p — ie (p — q't) 2 — ie 
x2 [L 2 (e, Z) (f + 75 ) + R 2 (e, Z) (f - 75 )1 u(ki) . (C.17) 



In the on-shell renormalization scheme, the t-channel counterterm part in 
equation (C.16) originates from the counterterms of two vertex and one elec- 
tron self energy loops. It reads 



B 



a/3 



Cf z = -Jf hp + 4P IR - 4) 

Qt 



(C.18) 
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with B^y from equation (C.10) and the ultraviolet and infrared poles 



P = + - Ie + In— - ln(2v^F 

n — 4 2 jj, v 



pIR 



// 1 2 
1 1 



n-4 2 



+ - 7 £ + ln— -ln(2v^r) 



/i=m e 



/i=m e 



n = 4-e, (C.19) 



n = 4 + 5. (C.20) 



Here, 7^ is Euler's constant. 

The w-channel virtual QED matrix element M^ z is obtained from Ai^z by 
the interchanges 



vi <-> u 2 , 



34,Z- 



(C.21) 



The symmetry (C.21) implies that, after the angular integrations, M.^ z * s 
not explicitly needed any more, because the terms corresponding to the in- 
terference of Ai^z with the tree level matrix element is obtained from the 
interference of M^,z with the tree level matrix element by interchanging s\ 
and S2- 



The virtual initial state QED matrix element for the CC3 t-channel is: 

M V _ 9PP' 9aa> yaf3 yrff M c 



D w (si) D w (s 2 ) 



't,W Ivl 12,W 1V1 U,W 



(C.22) 



with 



V af3 - e 1 
v t,w — e 



yap 
v t,aux 



1 n 



(4-n) 



(2tt)™ 



yaf3 
v t,box 



H=m e 



(C.23) 



where 

ya{3 



/-tap j n 
yap _ ^t,W _ ■ (4-n) / a P (yap , yap , yap N 
^t.aux - 1 ^ J \ V t,vertl + ^,vert2 + "t,self, 



. (C.24) 



H=m e 



The virtual and counter term contributions are those of (C.17) to (C.20) with 
corresponding Z couplings replaced by W couplings. The electron mass which, 
in principle, occurs in the NC2 t-channel propagators is absent in the charged 
current case where a neutrino is exchanged in the t-channel instead. Further, 
the w-channel diagrams are absent. 
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The virtual initial state QED matrix element for the s-channel annihilation 
represented by the amputated Feynman diagram in figure 4 is given by 



M 



v 



077' 0/3/3' 9aa> T/7 ' 



D v o(s) D w (si) D w (s 2 ) 



M^ w M& w g 3 (V°) T^(C.25) 



with 



T/7 

V s,V° 



2 I O s,V 

! ' ' < 16^ ~ 1 " 



i / d> ^ 7 



(27r) n 



s,vert 



fj.=m e _ 



(C.26) 



and 



T/7' 



s,vert 



= m(-^ 2 ) 



-2^- 7 ^)7 7 ' (2^-^ 7") 



p 2 + 2/c 2 p — ie p 2 — 2kip — ie 



p z — \e 



x 



L(e,\/ )(l+7 5 ) + ^(e,\/ )(l-75) u(h). (C.27) 



The s-channel counterterm part is a mere vertex counterterm given by 



Cl'yo = B]' y0 (2P + 4P IR - 4) 



(C.28) 



with B]y from equation (C.9). 



C.3 Initial State Bremsstrahlung Amplitudes 



As for the case of virtual initial state QED corrections, initial state brems- 
strahlung matrix elements are obtained by altering the amplitudes in equa- 
tions (C.3) to (C.7). From the Feynman diagrams of figure 7, one deduces the 
t-channel Bremsstrahlung matrix element for the NC2 case: 



M 



R,\ 
t,Z 



9/313' 



9o 



Dz(si) D z (s 2 ) 



'TrM^M^eRf/e^p) 



(C.29) 



with 



u(-ko 



Z2 
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x2 



L 2 (e, Z) (1 + 75 ) + R 2 (e, Z) (1 - 7s ) u(fci) . (C.30) 



An analogous result is obtained for the NC2 w-channel initial state bremsstrah- 
lung matrix element, 



M 



R,X 9aa> 9 '0/3' 



u,Z 



r>a/3n 
n u,Z 



D z {si) D z (s 2 ) 12 > z u ' z u ' z 
u{-k 2 ) < 7 7 — - - — 



u 2 



u 2 Ui 

z 2 Ui j 

x2 [L 2 (e, Z) (1 + 75 ) + R 2 (e, Z) (1 - 75 )1 u(h). (C.31) 



For the s-channel bremsstrahlung matrix element of the CC3 process one has 
the Feynman diagrams of figure 5 and obtains 



Ms > v ° - D v o(s) D w ( Sl ) D w (s 2 ) Ml2 > wMu > w93[V ]1 eR sV° £ ^ 

(C.32) 



with 



R% = u{-k 2 ) { ^ 



, 2k{-i^ 7^ -2^ , 



+ 



-7' 



z \ z 2 
x [L(e, V°) (1 + 75) + R(e, V°) (1 - 75)] u(h). (C.33) 



All the initial state currents introduced so far fulfill the condition of U(l) 
current conservation, 



p.RTvo = p,Rf" = p,Rf" = 0. 



(C.34) 
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There is no CC3 w-channel diagram. For the CC3 t-channel diagram, the sit- 
uation is slightly more complicated. Naively, only the two left diagrams in 
figure 7 contribute to the CC3 t-channel ISR corrections. For these two, cur- 
rent conservation is violated and may be established by adding as an auxiliary 
current the additional Feynman diagram with radiation from the t-channel 
propagator [7] (cf. the discussion in section 3.2). Then, the charged current 
t-channel matrix element is just given by (C.29) with Z couplings replaced by 
W couplings. 

In the above matrix elements the radiatively changed Mandelstam variables 
were introduced. In terms of the five-particle phase space parametrization from 
appendix B they are given by 

a\ + b\ cos6 R + c sin^R cos<f>R, 
a\ + b 2 cos9 R + c sin#R cos0#, 
a" — bi cos9 R — c sin#R cos<fi R: 
0% — b 2 cosOr — c sin#R cos0r. (C.35) 

Here m v = m e for NC2 and m v = m v = for CC3. In (C.35), the kinematical 
parameters 



a\ = ml 


- m 2 e - 


Sl 


+ 


s' + S , 

As' y- 


*+ - sLp cose) , 


a\ = ml 


- m 2 e - 


S 2 


+ 


s'-5 , 

As' y 


'+ + si (3 cos6>) , 


a\ = m 2 v 


- m 2 e - 


S 2 


+ 


s'-6. 
As' V 


/ + - s'_(3 cos6) , 


a-2 = ml 


- m 2 e - 


Sl 


+ 


s' + 5 , 
As' V 


t' + + si (3 cose) , 



h = ^ (s> + f3 cose -s'_), 

fx 



(3 sine (C.36) 

are derived from the representations (B.4) and (B.10) in appendix B. The 
variables 5, si, si, viV are introduced in sections 2 and 3. 




c = 
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The two invariants related to radiation from external legs are 



z\ = —2 kip 



z 2 = -2k 2 p 




(C.37) 



The equations (C.36) and (C.37) are exact. In (C.37), one has to keep the 
electron mass, since it gives rise to mass singularities. The electron mass in 
(C.36) may be neglected, but it was retained in some integrals for the sake of 
numerical stability. 



D Phase Space Integrals 

In this appendix, we collect the analytical solutions of the integrals needed for 
the semi-analytical cross-section calculation of processes (1.1). These integrals 
are necessary to obtain the twofold differential cross-section d 2 a/ dsi ds 2 in the 
case of virtual ISR and the threefold differential cross-section d 3 <7 / ds' dsi ds 2 
for initial state bremsstrahlung. Integrals are presented in the ultrarelativistic 
approximation, i.e. the electron mass is neglected wherever possible. The inte- 
grals have been numerically checked. The azimuth around the beam yields a 
factor 2tt. We do not elaborate on the pure s-channel in this appendix, because 
it is well-known that its initial state radiation is fully covered by the universal 
ISR corrections. 



D. 1 Fermion Decay Angle Tensor Integrals 

Integrations over final state fermion decay angles are carried out with the 
help of invariant tensor integration. Using the decay matrix element M^ 2 V 
from equation (C.12) one obtains 



Spins 



dr 12 Mf 2y Mr 2 

,v 




Spins _i q 



yf\(si,ml,m%) 4tt 



L^f^V) + R\f u V) 



8si 3 



x Sl ff 



+ K+^) (P1+P2) 



(D.l) 
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with L(f\, V) and R(fi, V) from (A. 2) and (A. 8). An analogous result is valid 
for the integration of M^ iV M^ iV over dr 34 . The result of the above decay 
phase space integral factorizes in the cross-section. Therefore one can use 
formula (D.l) to integrate the final state fermion decay angles for both the 
2^4 and the initial state bremsstrahlung 2— >5 particle phase space. 



D.2 Loop Integrals 



For the computation of the virtual initial state QED corrections to the studied 
processes, apart from the integration over final state fermion decay angles, two 
sets of integrals are needed. The first set is presented in this subsection and 
consists of the integration over the loop momentum p: 



with |n-4| < 1. (D.2) 



fi=m e 



The second set, the integration over cos-# (in 2 4 kinematics), is given in 
section D.3. It should be noticed that we present only the loop integrals for 
the virtual s- and fr-channel graphs, because, in the twofold differential cross- 
section, the interferences of the virtual w-channel graphs with the tree level fr- 
aud w-channel graphs equal the interferences of the virtual fr-channel graphs 
with the tree level w-channel and fr-channel graphs. This is easily verified by 
symmetry arguments. We use the definitions 
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p 2 — ie 
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(p — ki) 2 + ml — ie 


p 2 — 2pki — ie 
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(p + k 2 ) 2 + m 2 e — ie 


p 2 + 2pk 2 — ie 
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1 




n^ 


(p - qt) 2 - ie 





(D.3) 



and q 2 = t. For later use we define q 2 l = u. Dimensional regularization is used to 
consistently treat divergent loop integrals. The ultraviolet and infrared poles 
are taken from equations (C.19) and (C.20). For the sake of compact pre- 
sentation of the subsequent integrals we introduce two symmetry operations, 
namely 
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T 
U 



f(si,s 2 ) 



= /(s2,si) 
= g(-k 2 , -ki)- 



(D.4) 
(D.5) 



Further we introduce the notations 
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One obtains the following table of canonical integrals: 
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D.3 Virtual Corrections' Phase Space Integrals 



After the integration over the loop momentum, there is one non-trivial integra- 
tion left when the virtual initial state corrections to the NC2 or CC3 processes 
are evaluated. This is the integration over the boson scattering angle d in 
the center of mass system. We recall that the below integrals do not cover 
the interferences of virtual w-channel graphs with tree level t- and w-channel 
graphs, because these interferences are deduced from the interferences of the 
virtual t-channel graphs with the tree level w-channel and t-channel graphs by 
symmetry arguments. We introduce the notation 



dcos-# A 



'-■ma: 
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dt A 



(D.7) 



with t m i n and t max as defined in equation (3.4). 
Below, we use 



tl2 = s i + t 
Ua = s 2 + t 



(D.8) 



and adopt the notations of equations (3.4) and (3.6). The following list of 
integrals is obtained: 
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The above integrals 1) to 6) are familiar from the calculation of the tree level 
cross-section, while the integrals 7) to 29) genuinely originate from the virtual 
corrections. 
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D.4 Bremsstrahlung Integrals - First Series 



For the calculation of bremsstrahlung to processes (1.1), the phase space was 
parametrized as exposed in equation (B.3) and appendix B. After integration 
over the final state fermion decay angles, the first series of bremsstrahlung 
integrals is over the scattering azimuth and polar angles <pR and Or of the 
boson three-vector v\ in the two-boson rest frame: 



2tt +1 



[Mr = ^ l d ^ J dcos °R A ■ (D-9) 



We introduce a third symmetry operation, 



S 



f{cos8) = /(-cos0) (D.10) 



Taking into account that the contribution of the bremsstrahlung li-channel 
matrix element squared to the threefold differential cross-section equals the 
contribution of the bremsstrahlung t-channel matrix element, the following 
bremsstrahlung integrals are needed for the first series: 

1) [1} R = 1 2) [cos^r f(cos6 R )] R = 

3) [cos9 R ] R = 4) [cos 2 0/4 = V2 
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5) [cos 2 ^ = 1/3 



7) 



2s' 



R 



VX{s' + - s'_ cos9) 
2s' 



In 



_(s' +\/~>t)-2s' a+s'_8-{s' +8+\f>!)s'_ cos6>" 
_ («' - v 7 ^) -2s '<t+s'_ 8- (s> +8- s'_ cosB 



8) 
9) 
10) 

11) 
12) 
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d2 
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.) 
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R 



R 



si (s - 5 + - s'_ cos9) (s-S-VX-sL cos9) " MM 



ST 



RJ 



As' 

S2 d.2 d 2 
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15) 
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°12+ °12- 



= S 



ai_ + 6i_ COS6 1 = a 1+ + 6 1+ cosf? 

a 2 - + &2- cos# = a 2+ + 6 2 + cosf? 

a-s5 2 + S -s'_ (s' + S- 2a) - VX (sS - S -s'_ 
-(VX + s> + 6) S -s'_ 

s 2 a - s5 2 + ^s'_ (s' + 5- 2a) + \^X (s5 - ^s'_ 

+ (yx -s'-s) y_ 



s 2 a - s5 2 + (s'-S- 2a) - ^/X (sS + U_ 



= -(yx-s'+6) y_ 



= s 



V - s5 2 + -sL (s'-5- 2a) + VX(sS + -s'_ 
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52 



16) 
17) 
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At 

Bi 



Cl2 

B{ 



TDU 
D 2 



VC^ b d (a{ -h) + C 12 + B\ (a sl + b d ) /2 
VC^2 b d (a« - b 2 ) + C 12 + B% (a sl - b d ) /2 
^ &d ( a « + &2 ) + Cl2 + £ 2 " (a sl + b d ) 12 

(^) 2 K co^ + 5 12 ) 
-s'_ 2 Si (s - Si) 

s'} s (s — Si) — 2ssi (s — a) (s — s x ) + s 2 (s — a)' 
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-2 ffe 2 (a<& 2 + a^6i) + a si c 2 
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20) 



Mui. 



R 



T 



-tlU 2 . 
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29) 
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D.5 Bremsstrahlung Integrals - Second Series 



Below, the integrals needed for the last analytical integration in the brems- 
strahlung case are given. This integration is over the photon scattering angle 
6. Using the ultrarelativistic approximation and the notation 



[A\ e = -J dcos6 A, 

-1 



(D.11) 



the following integrals are required for the "second series" of bremsstrahlung 
integrals: 



1) [1]« 



1 
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1 + 



s'-cf 



2s! 



■ A' 

L c i + Try- -^c3 + L c5 
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In 1 + 
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5 
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s'(a - 5) 



s'(a - 5) 



13) 
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VZi J 
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— L c5 — L 
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2s(sL+5) + s'(a-5) / I 2s(sL+5) + s'(a-5) 
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A comment on the treatment of arguments of logarithms and dilogarithms is in 
order here. The logarithm is undefined for negative real values, the dilogarithm has 
a cut for real numbers larger than 1. This means that, in principle, infinitesimal 
imaginary parts from the propagator denominators have to be carried along in the 
whole calculation. This, however, is an unnecessary nuisance, because all integrands 
are real and regular. If infinitesimal imaginary parts were needed, because logarithm 
or dilogarithm arguments lay on cuts in final results, they were attributed to the 
invariant masses s\ and S2- As integrands are real and regular, this is a correct 
treatment, because it is then irrelevant for the integrand how an infinitesimal imag- 
inary part is entered. This technique will be used in subsequent integrals without 
further notice. 
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8 Integrals 36) to 41) were not directly calculated from the result of the Qr in- 
tegration. Instead, a Feynman parametrization was used to "linearize" l/£i« 2 and 
\jt 2 u\: l/t\u 2 = Jq da/ [t\a + u 2 (l — a)] 2 . Then, Qr, 9, and finally the Feynman 
parameter a were integrated. 
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In the calculation of the integrals the following quantities appear: 

A 3 = si (4 - 2a) 2 - /3 2 sA' 
£3 = si (s| - 2a) (s - 5) 

|ss' — ss 2 — s'si — 2si (s — it)| 



A 



2^3/4 = /5 



-4m 2 { Sl (s-(5) 2 + A' (s-si)} 
-S3 ± 



In principle, depending on the values of x 3 and X4, four cases must be 
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distinguished: 



(i) x 3 / 4 G M; |x 4 | > \x 3 \ > (3 

(ii) X3/4 G ]R; \x 4 \ > (3 > \x s \ 

(Hi) X3/4 G M; (5 > |x 4 | > |x 3 | 

(iv) A 3 < 0; x 3 = x% G C 



However, since the integrand is regular, it is sufficient to present the solu- 
tion for case (i). The solutions for the other three cases are then obtained 
by analytical continuation. In case (iv), which is relevant for only a very 
small fraction of the phase space, the problem of a numerically correct 
treatment arises. This problem was solved by computing the integral's 
value for a very nearby phase space point so that the expression of case 
(i) could be used. Having in mind that the integral is regular, it is clear 
that thus only a negligible error is introduced. The expressions X , X ± , 
and X 2 for case (i) are presented below. 
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